Field theory of the quantum kicked rotor 
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The quantum kicked rotor (QKR) is investigated by field theoretical methods. It is shown that the 
effective theory describing the long wave length physics of the system is precisely the supersymmetric 
nonlinear cr-model for quasi one-dimensional metallic wires. This proves that the analogy between 
chaotic systems with dynamical localization and disordered metals can indeed be exact. The role of 
symmetries is discussed. 
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Quantum mechanics tends to suppress the chaoticity 
of classical dynamical systems. The investigation of this 
phenomenon in periodically driven systems, i.e. systems 
that are governed by a Hamiltonian with periodic time 
dependence, has led to the discovery of one of the most in- 
triguing parallels between the fields of nonlinear dynam- 
ics and disordered solids: In the quantum kicked rotor 
(QKR), a typical representative of this class of system, 
the quantum mechanical suppression of chaos exhibits 
striking similarities to the phenomenon of Anderson lo- 
calization in disordered metallic wires 0. 

By definition, the kicked rotor is a point particle that 
moves freely on a circle. The particle is kicked period- 
ically in time, where the kick strength depends on the 
angular position. When a kick strength parameter fcr 
exceeds a certain threshold value, the dynamics becomes 
globally chaotic. In a statistical physicist's language, the 
chaoticity of the motion manifests itself as follows: An 
ensemble of particles prepared at time t = so as to 
have definite angular momentum Iq but arbitrary angu- 
lar coordinate 9, will diffuse in Z-space around the ini- 
tial condition Iq. In the corresponding quantum system 
(6 —> 0, I — ► /, [1,6] — —ill), the unbound diffusion 
in l-sp&ce is suppressed by localization. Numerical Jl]] 
and analytical |^,|| studies have shown that the QKR- 
localization is analogous to the Anderson localization dis- 
played by metallic wires with many channels (quasi Id 
wires). In particular, it has been demonstrated numeri- 
cally j| that a phenomenological modeling of the QKR 
by random band matrices can successfully explain essen- 
tial features of the localization phenomenon. Random 
matrix models of the same type are known ||, in turn, 
to describe the universal large distance physics of disor- 
dered wires. 

However, the equivalence between the rotor and quasi 
Id wires still has the status of a conjecture. A rigor- 
ous answer to the question whether this analogy is com- 
plete (and not just restricted to the bilateral appearance 
of localization) has not yet been given. Can a simple 
one-dimensional driven system indeed exactly mimic the 
behavior of disordered electronic conductors, which in- 
cludes a variety of complex phenomena that have recently 
been found |6|,^2|||? In the present Letter we are going 



to show that the answer is positive. This is done by 
mapping the kicked rotor onto the very same supersym- 
metric nonlinear cr-model that is known to describe the 
long wave length physics of disordered wires. The fact 
that both models can be described by the same effective 
field theory implies that all that is known about the QKR 
applies to disordered wires and vice versa. 

The QKR is defined by the time dependent Hamilto- 
nian 

p. A 00 

H = — + kcos(9 + a) 5(nr-t), 

n— — oo 

where the particle's moment of inertia has been set to 
unity and a G 1Z is a symmetry breaking parameter whose 
meaning will be explained below. To elucidate the anal- 
ogy between this Hamiltonian and disordered electron 
systems, one may consider the discrete time analog of a 
four-point Green function in (angular) momentum space: 



(h\G + (uj + )\l 2 )(l 3 \G-(^)\k) 



(1) 



±11-1 



where G ± {lu±) := Y%=o U n e luJ±nT = [l-(C/e M±T ) 
U = exp(zi 2 r/4) exp(ik cos(9 + a)) exp(iZ 2 r/4) denotes 
the Floquet operator, i.e. the unitary operator govern- 
ing the time evolution during one elementary time step, 

u>± = u a ± (u/2 + iO), and (. . .) Uo := r J^ /t dw {. . .)/2n 
is an average over the rotor's quasi-energy spectrum. 

Before turning to the quantitative discussion of the sys- 
tem, let us explain the meaning of the symmetry breaking 
parameter a. The localization length of a disordered wire 
depends on the behavior of the Hamiltonian under the 
time reversal transformation T: t — > —t, p — > —p, x — > x, 
where p and x are the momentum and the position. In 
the case of the rotor, where localization takes place in 
momentum rather than coordinate space, this symmetry 
operation is irrelevant. However, it has been shown || 
that the transformation T c : t — » —t, 6 —> —8, I — > I plays 
a role analogous to T in disordered metals. (Note that 
T c differs from T just by the exchange of momentum and 
position.) To couple the system to a T c breaking pertur- 
bation in a simple way, we put it on an angular lattice of 
spacing 1/(2%L), L G M, thereby giving it the topology 
of a ring of circumference L in momentum space JTG] . It 
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will turn out that the symmetry breaking parameter a 
then acts like a T c breaking Aharonov-Bohm flux pierc- 
ing the ring. 

In the following we deal with the correlator ([!]) by field 
theoretical methods. The strategy of our approach is dic- 
tated by the experience gained from both the analysis of 
disordered metals JO] and a recent field theoretical ap- 
proach iQ to Hamiltonian chaotic systems. Owing to 
the different formulation of periodically driven systems, 
however, the actual computational scheme deviates sig- 
nificantly from these cases. To simplify the notation, we 
temporarily focus on the case of unbroken T c symmetry, 
a = 0. 

Invariance under the transformation T c , which acts as 
an anti-unitary operator in the quantum system, results 
in the Floquet operator being a symmetric matrix when 
represented in the Z-basis (from now on we refer to all 
operators in Z -representation) . This makes it possible 
to decompose U = {(l\U\l')} by (Ue*" ±T ) ±1 = V±V£, 
where V± does not possess any symmetries other than 
unitarity. We choose V± = e ±iu±T/2 K%r£ /2 , where 

K a ±k := {(Z|exp(±if 2 T/4)exp(±ifccos(6> + a))\l')}, and 
write the Green functions appearing in (0) as jl3| 

G ± ^±) = ( yT ^ ) =: G±(w±) u . (2) 

In the next step we introduce a superfield ip — {ip\ a u}, 
X,a,t — 1,2, with complex commuting (anticommut- 
ing) components 4>a=i{^a=2), and consider the gener- 
ating functional 

yp(^,^)exp[-^(G- 1 + J)V] , (3) 

where G = £ AR <g> 1 BF <g> G+ + Ef K ® 1 BF (g> G~ ) , 
matrices with subscript 'AR' ('BF','T') act in the two- 
dimensional spaces of A (a,t) indices (the ^-indices refer 
to the matrix structure appearing in (|2|)) and (E^ := 
8u'8jji , X = AR, BF, T. Here and below, indices that are 
not indicated explicitly are summed over. Expressions 
like (0) can readily be obtained from (g) by differentiat- 
ing twice with respect to matrix elements of the source 
field J. As we are interested in the general structure of 
the theory, rather than in the calculation of any particu- 
lar correlation function, we henceforth omit J. 

After a few elementary manipulations, namely matrix 
transpositions and regrouping of integration variables, 
the Gaussian integral (^) takes the simple form 

J V{(j),x)e~^ + * x)+ ^ El ^® V+x+ * E ^® V - 4 ' ', (4) 

where the fields <f> — {4>\ a tl} an d X — {x*ati} comprise 
components of both ip and tp. Instead of displaying the 
structure of these new quantities explicitly, we merely 
note two essential features that fix their functionality as 



integration variables: (i) 4> and x are independent of each 
other and (ii) they possess the symmetry Y = (Iar <8> 
M)Y, Y = <f>,x, where M = E^ F ® a\ + E^ F ® (icr|) 
and a l x (i = 1,2,3, X = AR, BF, T) denotes the Pauli 
matrices. 

The next step in the construction of the field theory is 
the average over the phase exp(jwoT), which plays a role 
similar to the energy average employed in Ref. In 
that case, energy averaging led to a quartic (~ (tpifj) 2 ) 
non-local contribution to the action of the field theory. 
The latter was eliminated by means of a matrix valued 
auxiliary field Q that coupled to the dyadic product ipip. 
The phase average to be carried out in the present prob- 
lem produces in addition to the quartic term an infinite 
series of higher contributions to the action. We have 
not succeeded in decoupling these terms by elementary 
means. On the other hand, the experience gained from 
previous diagrammatic analyses || of the QKR suggests 
that a field coupling to ip-ip should again describe the large 
scale physics. 

The problem of identifying this field is solved by a re- 
cently discovered identity O] that adapts the Hubbard- 
Stratonovich transformation to averages over unitary op- 
erators. In the special case under consideration, namely 
a phase or U(l) average, this identity reads: 

(efavi+v***^ = J X> M (Z,Z)e^i z ^+'l2^i ; (5) 

where u := exp(iu T), m = V+| Wo=0 Xi, m = X2Vl\^ ( _ , 
all subscripts refer to the A-indices ('AR'-space), Z = 
{Z atl, a? vv} is a non-local (in I) 4 x 4 supermatrix field, 
T>n{Z, Z) = V(Z, Z)sdct(l-ZZ) with 'sdet' the superde- 
terminant, and J T>(Z, Z) stands for the integral over the 
matrix elements of Z and Z :— Z T cr B p. 

The proof @| of (|) makes use of group theoretical 
concepts and the theory of generalized coherent states 
and is too lengthy to be reported here. We note how- 
ever that the field Zip takes values in the unrestricted 
set of 4 x 4 complex supermatrices. The latter can be 
interpreted as a space parameterizing the coset space 
G/K, K = {k G G\ka\ K = a\ K k) C G, where G is 
the group of 8 x 8 supermatrices g subject to the con- 
straint g^rjg = ry, 77 = (er AR ® E^ F + 1 AR ® Eg F ) ® 1 T . 
This coset space is the field manifold of a 'unitary' su- 
persymmetric cr-model that is twice as large as in the 
usual case 0] on account of the extra T-space indices. 
The relationship between the Z-field and this manifold 
is the first indication of the fact that we will end up with 
a nonlinear cr-model. 

After this comment on the formalism, we proceed to 
apply (||) to the construction of the field theory for the 
rotor. To that end we insert (g) into (^) and perform the 
Gaussian integration over the fields <f> and %■ As a result 
we obtain for the generating functional (at J = 0), 
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J V(Z, Z) cxp str 



ln(l - ZZ) - ±ln(l - Zt- x Z t t) 
-±ln(l-e iUT ZU T- 1 Z T TU%) , (6) 



where Uq = U\ a _ Q , r = M<7g F and the supertrace 'str' 
includes a trace over the Z-space. So far all manipula- 
tions have been exact. We next restrict the field theory 
to its infrared limit, which describes the long time/large 
distance physics we are interested in. 

The action of the field theory (Q) vanishes for fields 
Z proportional to unity in angular momentum space if 
lo — > and the constraint 



Z = r~ 1 Z T T 



(7) 



is imposed. Field configurations that violate this sym- 
metry are 'massive' and cannot contribute to the long 
range correlations of the model. We therefore restrict 
the integration in (||) to the field manifold specified by 
(0). (The integration over the massive quadratic fluc- 
tuations around this manifold yields a factor of unity 
by the supersymmetry of the model.) Note that the 
unitary coset space G/K subject to the constraint ([?]) 
defines the field space of the ^orthogonal nonlinear a— 
model. As a preliminary result, we thus find that our 
field theory has the same symmetries as the one de- 
scribing time reversal invariant disordered metals. What 
happens when this symmetry is gradually broken by 
the introduction of a small finite value of a? In that 
case, the decomposition of the time evolution operator 
has to be generalized to (Ue iuJ±T ) ±1 = V±V± aT , where 
ya = e ±iu±r/2 K a Au lurtrie r steps can b e rep eated 
in essentially the same way as before and we again ar- 
rive at (||). The only change is that the fields appearing 
in the action have undergone a 'gauge transformation' 
Zi\i — > Z%, := exp(-iala^) Zui exp(zaZ'cr^), and similarly 
for Z. 

To carry the analogy to disordered metals further, we 
need to expand the action around the limit Zu> = SwZq, 
lo = 0. Details of this somewhat tedious calculation 
will be presented elsewhere |l6} |. Here we restrict our- 
selves to a rough sketch of the main ideas. We first 
subject the action in (^J) to a 'semiclassical approxima- 
tion'. The expansion parameter of this approximation 
is h/Sl <C I, where SI is the typical angular momentum 
scale over which the relevant Z-fields fluctuate jlT]]. As 
a result, (i) the symbol 'str' in (^J) no longer includes a 
trace over angular momentum space but rather an inte- 
gral over the phase space coordinates (1,8) of the clas- 
sical rotor, (ii) Z(l,l') is replaced by its Wigner trans- 
form Z(l,0) (we temporarily suppress the superscript 
a in Z a ), and (iii) U ZUq is replaced by Z u , where 
f u (l,8) = f(9 + r(l + k sin 9),l + k sin 9) denotes the clas- 
sical one time step evolution (standard map) of a phase 
space function. 



The angular variable 9 of the standard map is a rapidly 
relaxing degree of freedom, which leads us to expect that 
only ^-independent field configurations contribute to the 
long wave limit of the model. To formulate this state- 
ment in a quantitative manner we do a Fourier transform, 
Z(l, 9) — J2m=-<x> Zrn(l) exp(im9), and observe that the 
non-zero modes Z m ^o are 'massive'. Integration over 
these fields in Gaussian approximation yields an effective 
field theory for the massless zero mode Z(l) := Zq(1). 

We finally expand the action S[Z] in terms of slowly 
fluctuating fields. This program is carried out most eco- 
nomically in the Z-Fourier space, i.e. in angular coor- 
dinates fyl. The small parameters of this expansion 
scheme are lot and the characteristic 'momentum' <fi (<fr 
is an angular variable) of slowly fluctuating fields Z(<j>), 
which is of order (j> <C k^ 1 <C 1. Concerning the former, 
we note that the mean quasi-energy level spacing of the 
model is A = 2tt(Lt)~ 1 . Since we are interested in small 
frequencies of 0(A), we have lot ~ L _1 <C 1. To leading 
order in lot and <f> the action of the generating functional 
reads 

S[Z a ,Z a ]~ j dx[- l -^str(l-Z a (x)Z a (x))- x - 
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d x ,d xS tT\n(l-Z a (x)Z a (x'))\ x , =x , (8) 



where we have Fourier-transformed back to angular mo- 
mentum space, taken a continuum limit (i.e. the variable 
x is a smoothed version of the Z-index, ^ — > J dx) and 
D = k 2 /2 + ... is the classical diffusion coefficient of the 
rotor jl9| . (The dots indicate oscillatory corrections pfj] 
to D that result from elimination of the non-zero modes 
and are smaller than the leading term by powers of k.) 
Introducing an 8 x 8 matrix field Q by 
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we can rewrite the functional integral as 



VQexp J str (J^V a QV a Q + ^Q4r 
V a = V + ta[«4. ■ ]. 



(9) 



which is precisely the nonlinear tr-model for a quasi one- 
dimensional metallic ring in the presence of a T-breaking 
(Aharonov-Bohm) vector potential of strength a. 

Because of its importance for an understanding of the 
localization physics of wires, the model (H) has been in- 
vestigated thoroughly Let us now review some of 
its essential properties in the terminology of the kicked 
rotor. For times less than 0(Tk 2 ) the kicked particle per- 
forms a diffusive motion in momentum space. On larger 
time scales, quantum localization confines the particle to 
stay within a volume specified by the localization length 
£ = fc 2 /2. For 'flux' strengths a ~ L^ 1 , the orthogo- 
nal symmetry of the model is broken and one expects a 
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doubling of the localization length £ D — > = 2£ Q [ p"l| |. 
(Note, however, that a max = 2-7r/L corresponds to one 
'flux quantum' O penetrating the ring. As the physics 
of Aharonov-Bohm geometries is c/>o-periodic, a max is 
the maximum field strength that can be realized in our 
model. For systems with L > £ Q this strength does not 
suffice to cause a crossover from £ Q to 

A careful look reveals that the localization length £ Q 
predicted by our analysis is four times larger than the 
length found in numerical work (cf. e.g. Ref. 
We believe that this discrepancy is caused by an am- 
biguity in the convention of what is called a localiza- 
tion length: The 'field theoretical' localization length 
determines the exponential decay of the average tran- 
sition probability (|G + (w + ; I, ~ exp(— \l — l'\/£ ) 
between two remote states \l — l'\ ^> £ G . In numerical 
measurements, however, one computes the average of in- 
dividual decay constants, which is to say that one calcu- 
lates \l - l'\/£ n = (-hi|G+(w + ;Z,Z')| 2 ) wo . The lengths 
£o and thus defined do not coincide in general. Ac- 
cording to Ref. |^] they are related by £ = 4£„ for 
quasi lei-wires. In view of this, our analytical result does 
agree with the numerics. To summarize, we have mapped 
both the unitary and the orthogonal quantum kicked ro- 
tor on the supersymmetric nonlinear cr-model for quasi 
Id wires. This proves the longstanding conjecture that 
the universal properties of these two classes of system 
are indeed the same. Our mapping is straightforward 
and direct and avoids some approximations made in ear- 
lier work, namely (i) the replacement of the deterministic 
rotor by a stochastic model and (ii) the passage from uni- 
tary to Hermitian randomness. Note that the rotor-metal 
analogy is not restricted to the phenomenon of strong lo- 
calization. It has recently been shown that quantum in- 
terference in metals manifests itself in various pre-stages 
of localization such as non-trivial wavefunction statis- 
tics or t ne appearence of pre-localized states [||. 
The exact correspondence between quasi Id wires and 
the rotor suggests that these effects must be observable 
in the latter, too. In fact, the rotor may be an ideal 
model system for highly accurate numerical analyses of 
these pre-localization phenomena, since it can be imple- 
mented more efficiently on a computer than can weakly 
disordered multi-channel wires. 

We have benefitted from discussions with O. Agam, S. 
Fishman, P. Frcche, M. Janssen and A.Mirlin. 
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